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Abstract 

We propose an efficient mechanism for the operation of writing spin in a 
quantum dot, which is an ideal candidate for qubit. The idea is based on 
the Andreev reflection induced spin polarization (ARISP) in a ferromagnetic 
/ quantum-dot / superconductor system. We find that on the resonance of 
Andreev reflection, the spin polarization of quantum dot strongly denpends 
on the magnetization of ferromagnetic electrode, and the sign of the spin po- 
larization is controllable by bias voltage. In the presence of intradot Coulomb 
interaction, we show that ARISP effect can still survive as long as the charging 
energy U is comparable to the superconducting gap A. Detailed conditions 
and properties of ARISP are also discussed. 
PACS numbers: 74.50.+r, 73.63.Kv, 72.25.Dc, 85.35. Gv 
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Qubit is the basic unit in the achievement of quantum computing. Among many quantum 
two-level systems, the spin of quantum dot (QD) is an ideal candidate for such purpose [1]. 
Firstly, QD fabricated in semiconductor 2DEG can be well controlled by metallic gates, 
parameters such as resonant levels, tunnel barriers, and charging energy are experimentally 
tunable. Second, the spin coherent time is extremely long in semiconductors (exceeding 
100ns at low temperature [2]), which is orders larger than charge coherent time. Moreover, 
due to 0D nature of QD, many spin-flip mechanisms are further suppressed , resulting in an 
even longer lifetime of spin in QD [3]. 

One of the challenge to exploit the spin of QD as qubit is to efficiently control the local 
spin polarization in QD. A natural idea is to apply a magnetic field and induce the spin 
polarization by Zeeman splitting [4]. But the scheme contains practical difficulties because 
the field is required to be of the order of Tesla but confined within the scale of QD. An 
alternative idea is to attach QD to ferromagnetic electrodes [4,5]. In principle, the asym- 
metry between two spin bands in the electrodes may induce spin-dependent renormalization 
and broadening to the resonant level of QD and lead to spin polarization. However, such 
mechanism requires relatively large magnetization in electrodes and relatively strong cou- 
pling between QD and electrodes. More important, the sign of spin polarization in QD can 
not be controlled electrically. 

It is the purpose of this paper to present an efficient mechanism for writing spin in a QD 
with ferromagnetic (F) and superconducting (S) electrodes. The idea is motivated by the 
fact that F can provide spin-polarized electrons while S can only accept Cooper pair which 
is spin singlet. We shall show below that on the resonance of Andreev reflection (AR) [6-9], 
the spin polarization of QD strongly depends on the magnetization of F, and the sign of spin 
polarization is controllable by the bias voltage. The effect can survive even in the presence 
of strong Coulomb interaction, as long as the charging energy of QD is comparable to the 
superconducting gap in S. 

non-interacting case-To see the physics clearly, we start with the non-interacting fer- 
romagnetic / quantum-dot / superconductor (F-QD-S) system modelled by the following 
Hamiltonian, H = H F + H s + H dot + H T , in which H F = J2ka( £ k — °~h)a) ka ak a is for the F 
electrode, H s = £ po . e v h\ja v<J + £ p (A&J T &t pl + H.c.) for the S electrode, H dot = Y.a E ^ <T c (y 
is for the noninteracting QD, and H T describes the tunnel couplings between QD and elec- 
trodes. We introduce the retarded, advanced, and lesser Green function in the Nambu 
representation: 
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The occupation numbers (nj) = (cjc-f) and = (cjq) are related to the integral of G < . 

In the non-interacting case considered here, G r , G a , and G< can be solved exactly in 
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Keldysh formalism as G r = (g„ 1 - ££j G a = (G r ) f , and G< = G r E<G a , where g is 
the Green function of the isolated QD (see Eq.(7)) and So the self energy caused by the 
tunnel coupling with electrodes. In the wide band limit, S can be evaluated analytically 
as S = Efo+Sso with 
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in which r^-f, and T s are the coupling strengths between QD and electrodes (notice 
that IVf 7^ F F i), f(uj) = l/(e w / fcsT + 1) is the Fermi function, is the bias voltage between 
F and S [10], p r (id) and p < (uo) are the generalized BCS density of states defined as p r {oo) = 



Fig.l shows the curves of the occupation number (n a ) vs the resonant level E = = 
E±, with the ratio r^/T^f = |. (a), (b), and (c) are corresponding to the bias voltage 
V > 0, V = 0, and V < 0, respectively. For V — 0, there is a step from to 1 in 
the occupation curve, indicating an electron filling when E passes the chemical potential 
Pf = Ps- Moreover, the curves of (n-f) and are almost overlapped. For V > (V < 0), 
the step is shifted to £o = Pf = V", and there emerges a resonant dip (peak) pinned at 
Eo = Ps — 0. The results can be understood as follows (see also the schematic diagram in 
the bottom): Since no single particle states are available within the superconducting gap, 
electron filling to the QD is mainly determined by the F electrode, resulting in a step linked 
to the chemical potential of F. When the resonant level lines up with the chemical potential 
of S, however, two-particle process AR may occur, in which a spin| and a spinj electron in 
QD can leak into S by forming a Cooper pair and vice versa [11]. As a result, an Andreev 
dip (Andreev peak) is superposed on the step-like curve. The most remarkable features of 
AR dip (AR peak) are: (1) The spin polarization of QD, m = (n^) — strongly depends 
on the magnetization of F. (2) The sign of m is controllable by the bias voltage, m > 
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for V > and m < for V < 0. Hereafter, these nontrivial properties of AR resonance is 
referred to as Andreev reflection induced spin polarization (ARISP). Qualitatively, ARISP 
effect can be interpreted as follows: For V > and E = fi S , spinj and spinj electrons in QD 
tend to form Cooper pair and enter S. Since F can provide more spinj electron that spinj 
electron, the spinj electron will be depleted by spinj electron in the AR process, resulting 
in a strong spin polarization in QD. For V < and Eq = a Cooper pair is converted into 
a spinj and a spinj electron in QD. It is much easier for spinj electron than spinj electron 
to escape to the empty states of F, resulting in a reversed spin polarization. 

For quantitative analysis, we evaluate (n a ) near the AR resonance and obtain (nj) = 
1 — n(p) and (n^) = 1 — n(—p) for V > 0, (nj) = n(p) and (nj.) = n(—p) for V < 0, where 

° M 4"- p) 4Jig(i-^ + r r|(i- P 'l + r| - (4) 

in which p = (Tp^ — TpJ/lTp^ + ^fi) is the magnetization in F and Tp = (Tp^ + Tpi)/2 is 
the spin averaged coupling strength. Fig. 2a shows the AR dip for different magnetization p 
of F electrode (also true for AR peak if upside down the plot). Notice that the resonance has 
the Lorentzian line shape with half-width ^[T 2 F (1 — p 2 ) + T|] /4(1 — p 2 ), implying that AR 
resonance is broadened with the increase of p. Fig. 2b shows the maximum spin polarization 
m = m(E = 0) vs p for different coupling strength ratio r = Tp/T S - Contrary to the 
intuition, ARISP effect is most pronounced when the ratio r is small. Especially, in the 
limit r — > 0, (n-f)o = |(1 ± p) and (^j.)o = |(1 T V) (upper sign for V > and lower sign 
for V < 0). This means that both F and S electrodes are important in the ARISP effect: 
F provides the asymmetry between two spin categories, while S reinforces the asymmetry 
through AR process. Another noteworthy feature of ARISP is that the maximum spin 
polarization is determined by the ratio r and p, rather than the magnitudes of r^, F F j_, 
and Ts- Below we shall take into account the intradot Coulomb interaction to investigate 
whether ARISP can still survive. Due to the symmetry between electron and hole, only the 
case of V > is discussed. 

interacting case-To include the Coulomb interaction, we add the term Un^n^ to H dot , 
which makes F-QD-S a strong correlated system. Notice that in the limit U — > oo, double 
occupation is forbidden in QD, and electron number can not fluctuate by two. Hence AR is 
completely killed by Coulomb interaction in this limit. If, however, the charging energy U is 
comparable to the superconducting gap A, AR can still occur with the aid of bias voltage. 
Therefore, techniques for U — > oo limit (e.g., slave boson method) can not be applied to the 
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calculation of ARISP. Alternatively, we adopt the equation of motion (EOM) method, which 
is known to be reliable in the Coulomb blockade regime, and qualitatively correct for Kondo 
physics [12,13]. Moreover, EOM solution becomes exact in the U — > limit. In the EOM 
approach, one can derive the equation for the retarded Green function ({A(ti) \B(t 2 ))) r by 
differentiating with respect to t\ or t 2 , with new Green functions generated in the equation. 
To close the equation chain, we make the truncation in those Green functions containing 
two electrode operators in a mean field manner. After some algebra, the equation for G r 
can be obtained as 



AG r + G r A = 2N + B + B , 



(5) 



in which B = (g 1 1 —Y, r l )U 1 , A = B(g 1 — XJJ) + T, r 2 , A or B represents the transpose of A 
or B, U and N are defined as 
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Sq and H\ are the corresponding dressing self energies due to tunnel coupling with electrodes. 
£2 is the self energy contributed by the spin flip in the cotunneling process, which is related 
to Kondo physics. In the wide band limit, these self energies can be evaluated analytically 
[14]. Let £[= S^.+S^., (find T, r F0 and ££ in Eq.(2) and Eq.(3)) 
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in which Sl H ee -I J (f ), S2 (o;) = -I J (*) * S3 ^) = - J/ (I ) I , * 4 M = - JJ (f) 
andc^o = a>, a>i = uj — E^—Ei — U, uj 2 = w + -&r+-Ej. + t/, = w — £ , -j-+£'j_. The dimensionless 
functions /, J, Q are defined as 



arcsin x 



I(x) 



J(x) 



Q{x,y) 



X\f\ — X 

m 



7 ) 



2xVx^l 

7T 1 

2 V / W J 



|x|v? 



|x| < 1 

=5=1= In (\x\ + V^ 2 - l) , |x| > 1 
x\ < 1 



2 v^ 1 ! 

i / 1 
2 



stgn(x) 



+ 



\x\ > 1 



,e x + 1 e" + 1 

with ^ being the digamma function. 

As for G < , we invoke the stationary condition 
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where E 1 ^ = + £j_ + U and cr 2 is the 3rd Pauli matrix. We remove the integral over energy 
and write down the equation for G < as 



(S5 + \e u <j)j G< - G< (Eg + \e u o x 



= G r E<- 



-E<G a 



(17) 



which can be interpreted as "detailed" stationary condition, i.e., the observables in the 
energy interval uj and u> + du are time-invariant in the steady state. We note that Eq.(16) 
is exact while Eq.(17) is an approximation which neglects high order fluctuations in the 
presence of strong correlation (for details see [16]). Nevertheless, this approximation is 
sufficient for the study of ARISP which works in the Coulomb blockade regime. In addition, 
Eq.(17) becomes exact in the U — > limit. For comparison, we also try the Keldysh equation 
G < = G r E < G a , and employ the commonly used Ng's ansatz for E< [15,13]. Although the 
calculated results are qualitatively agree with each other, the convergence is much poorer 
in the latter approximation scheme. To sum up, Eq.(5) for G r and Eq.(17) for G< will be 
applied to the numerical study of ARISP. 
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Before presenting the numerical results, we qualitatively analyze the physics in the prob- 
lem. As seen in the non-interacting case, spin polarization occurs on the AR resonance. 
The conditions for AR resonance are: (1) QD is occupied with even number electrons and 
(2) a pair of electrons can transfer freely between QD and S. In the noninteracting case, 
these conditions amount to > I^r and E = hr for AR dip, or hl < jiR and E = jiR 
for AR peak. In the presence of Coulomb interaction, the occupation number of QD is 
quantized, and AR occurs in the even occupation region. Notice that the occupation of QD 
is mainly determined by F, let us cut off the S-QD coupling and consider the electron filling 
in F-QD system. QD has four occupation configurations: |0), ||), ||) and |||), with energy 
0, — fi F , Ei — fi F and + E^ + U — 2/ip, respectively. It is energy favorable that QD is 
empty when E — fi F > 0, singly occupied when —U < E — fi F < 0, and doubly occupied 
when E — fi F < —U. On the other hand, AR resonance requires energy degeneracy when 
moving two electrons from QD to S, so that Coulomb blockade is lifted. It is straight forward 
to write down the equality, E^ + E^ + U = 2fi S , meaning that a spinf and a spin j electron of 
QD enter S and form a Cooper pair at the chemical potential fis- Therefore the conditions 
for AR resonance are: E — fx$ — an d E — fx F < —U for AR dip, or E — fis = — ^ and 
Eo — Hf > for AR peak. The key point is that electron filling is linked to fip while AR 
resonance is linked to yU^, and fip ^ us in nonequilibrium. 

Fig. 3a shows the curve of (n a ) vs E for p = 0.5. As expected, a spin polarized AR 
dip is superposed on a step-like pattern. The steps from to \ and \ to 1 are located 
around E — fip = and E — fip = —U, while AR dip located around E — fx$ = 
which are in coincidence with the above discussion, (n-f) and (n^) are nearly equal off 
the AR resonance. This is because the weak F-QD coupling can only induce small spin- 
dependent renormalization and broadening to E<\ and E^. On the contrast, (n-f) and (n^) 
are dramatically different on the AR resonance, which relies on the ratio of IVj to Tp^ rather 
than their magnitudes. The up-right inset shows the detailed lineshape of AR dip for several 
p. Comparing with Fig.2a, one can see that Eq.(4) derived in the non-interacting case can 
also describe the ARISP in the presence of Coulomb interaction. In fact, Coulomb interaction 
plays the role of quantizing the electron number of QD. Coulomb blockade is lifted for the 
conditions discussed above, under which AR occurs as if through a noninteracting QD. Next, 
we discuss the choice of charging energy U, which is determined by the size of QD. Fig. 3b 
and Fig.3c shows the curves of (n a ) vs U on the AR resonance (E — [is = — f)- Fig. 3b is 
for the case of V < A, in which single particle process is forbidden due to superconducting 
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gap and AR process dominates. One can see in the plot that the spin polarization is almost 
independent on U when U < 2V, but gradually suppressed when U > 2V. The reason is 
as follows: on the AR resonance, QD favorites double occupation when U < 2V and single 
occupation when U > 2V. For double occupation, as shown in Eq.(4), ARISP is determined 
by the ratio p and r, and independent on other parameters. For single occupation, the 
effective S-QD coupling is greatly suppressed, and therefore ARISP vanishes. One tends 
to think that increasing the bias voltage V may help to overcome the charging energy U. 
This is true as long as V < A. When V > A, however, both single particle process and 
AR process are allowed, and the situation is more complicated. Fig. 3c shows (n a ) vs U for 
V > A. It turns out that the spin polarization is still independent on U when U < 2 A, 
but has a maximum at U — 2 A, suppressed when 2 A < U < 2V, and gradually vanishes 
when U > 2V. The suppression in the range of 2A < U < 2V can be attributed to the 
onset of single particle process. We note that EOM solution is valid for relatively small U, 
and the suppression of ARISP is probably underestimated when U > 2A. Nevertheless, the 
maximum at U = 2A is reliable, which is related to the singularity in the density of states 
of S electrode. 

In conclusion, we have proposed an efficient mechanism for writing spin in a QD which 
is based on the ARISP effect in F-QD-S system. The scheme has the advantages that 
the magnetization of F is not required to be strong and the sign of spin polarization can be 
controlled in fully electric manner. Calculation shows that the optimal conditions for ARISP 
are: T F < T s , U = 2A, E - /j s = -f , and \V\ > U/2. The properties of ARISP can be 
described by Eq.(4). In practice, the resonant level E can be tuned by gate voltage, and AR 
resonance can be monitored by the small tunnel current between F and S electrodes.. In the 
context of intensive research and impressive progress in F/2DEG, S/2DEG, and F/S hybrid 
structures, the proposed F-QD-S system should be feasible with up-to-date nano-technology. 
We are looking forward to hearing the relevant experimental response. 
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FIGURE CAPTIONS 



Fig. 1 The occupation number (n-f) (solid) and (raj.) (dotted) vs the resonant level E in a 
non-interacting F-QD-S system. The superconducting gap A = 1 is set as energy 
unit, the bias voltage V = 0.5, 0, and —0.5 for (a), (b), and (c), respectively. Other 
parameters are: T F = 0.01, T s = 0.1, k B T = 0.02, and p = 0.5. The diagram in the 
bottom schematically shows the AR resonance in non-interacting F-QD-S system. 

Fig. 2 Detailed analysis of the Andreev dip in Fig. la: (a) (n-f) (solid) and (dotted) vs 
E for different magnetization p of F electrode, (b) m = ((n-f) — (nj.)) Eo=0 vs p for 
different coupling strength ratio r. 

Fig. 3 (a) The occupation number (n-f) (solid) and (n^) (dotted) vs the resonant level E in 
an interacting F-QD-S system. The superconducting gap A = 1 is set as energy unit, 
the bias voltage V = 0.75, and the charging energy U = 1.2. Other parameters are 
the same as Fig.l. The up- right inset shows the details of the Andreev dip around 
Eq — — The middle inset illustrates the possibility of using ARISP effect to 
manipulate the spin states of a QD array, (b) and (c) show the curves of (n^) (solid) 
and (ni) (dotted) at E — us = — f vs the charging energy U, with the bias voltage 
V = 0.5 and V = 1.5, respectively. 
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